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The close analogy between cluster percolation and string proliferation in the context of critical phenomena is 
studied. Like clusters in percolation theory, closed strings, which can be either finite-temperature worldlines or 
topological line defects, are described by a distribution parametrized by only two exponents. On approaching 
the critical point, the string tension vanishes, and the loops proliferate thereby signalling the onset of Bose- 
Einstein condensation (in case of worldlines) or the disordering of the ordered state (in case of vortices). The 
ideal Bose gas with modified energy spectrum is used as a stepping stone to derive general expressions for the 
critical exponents in terms of the two exponents parameterizing the loop distribution near criticality. 



I. INTRODUCTION 

The quest for a geometrical description of phase transitions 
has a l0|Bg history going back to ideas first put forward by 
OnsageJll in the context of the A-transition in liquid ^He. The 
relevant geometrical objects in this transition are topological 
line defects. The description envisaged by Onsager is one en- 
tirely in terms of these 1 -dimensional objects, with their ge- 
ometrical properties such as fractal dimension and configura- 
tional entropy. The phase transition is characterized in this 
picture by a fundamental change in the typical string size. 
Whereas in the superfluid phase only finite strings are present, 
at the critical point infinite strings appear — similar to the sud- 
den appearance of a percolating cluster in percolation phe- 
nomena at criticality. 

The analogy between percolation of clusters andjarolifera- 
tion of strings has been noted by various authorsEoDu. As we 
will see, it derives from a similarity in the cluster and string 
distribution. Both have the same form containing two factors, 
one related to the entropy of a given cluster or string configu- 
ration, and the other related to the Boltzmann weight assigned 
to the configuration. Close to criticality, each of these factors 
is parametrized by a single exponent. One specifies the alge- 
braic behavior of the distribution at criticality, while the other 
describes how the Boltzmann factor tends to unity upon ap- 
proaching the critical point. Physically, the unity of the Boltz- 
mann factor implies that clusters or strings can grow and, thus, 
gain configurational entropy without energy cost. When the 
Boltzmann factor is not unity, the clusters and strings are ex- 
ponentially suppressed. In the context of strings, the transition 
between a phase consisting of finite strings only and one hav- 
ing infinite strings, is called a Hagedorn transition. We will 
refer to the appearance of infinite strings as proliferation. 

A central position in our arguments is taken by an ideal 
Bose gas with modified energy spectrum. The reason is that 
although noninteracting, the model has nontrivial critical ex- 
ponents, which are known exactly, while at the same time it 
can also be mapped onto a loop gas of worldlines in an ex- 
act way. The worldlines form closed loops because in the 
absence of external fields and sources they cannot terminate 
in the system. The map onto the loop gas allows us to con- 



nect the critical exponents of the phase transition to the two 
exponents parameterizing the worldline loop distribution near 
criticality. Using general scaling relations, these results can 
then be generalized to interacting loop gases representing sta- 
tistical models. Each universality class is defined by a loop 
distribution with specific values for the two exponents from 
which all the critical exponents follow. 

To close the circle, the resulting loop gas description of 
critical phenomena can also be applied to phase transitions 
involving (interacting) vortex loops — like the A-transition in 
liquid "'He. As for worldlines, vortices cannot terminate in the 
system and, therefore, form closed loops too. When the word 
string is used in the following, the two different 1 -dimensional 
objects worldlines and vortices should be kept in mind. 

As an aside, we believe that the theory of vertex loop gases 
discussed here has also a bearing on homogeneous superfluid 
turbulence. 

The paper is organized as follows. In the next section those 
essentials of percolation theory are recalled which later on in 
the paper become important to establish the connection with 
the proliferation of strings. In Sec. Q Bose-Einstein conden- 
sation in an ideal Bose gas with a modified energy spectrum 
is studied from the perspective of worldlines. In Sec. nu- 
merical work on random string networks is discussed in re- 
lation to uncorrected percolation, followed by a discussion 
of correlated percolation in Sec. In Sec. Q thermal phase 
transiti ons a re formulated in terms of proliferating worldlines. 
In Sec. VII, a similar formulation is discussed for phase tran- 
sitions involving pr olife rating vortices. The paper ends with 
conclusions in Sec. VIII. 



II. REVIEW OF SITE PERCOLATION 

In this sectio^, we briefly recall some basic aspects of per- 
colation theoryQ which are important for our purposes. Con- 
sider a lattice, with each lattice site being randomly occupied 
with a probability p, say, independent of its neighbors. Island 
of next-neighboring occupied sites form clusters. The prop- 
erties of these clusters as a function of p form the subject of 
percolation theory — or more precisely, of uncorrelated peico- 
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lation theory because the occupations of different lattice sites 
are uncorrelated. As p increases, clusters become bigger, and 
at some critical value pc, a cluster spanning the entire lattice 
first appears. Near this percolation threshold, various quan- 
tities show power-law behavior. As for thermal critical phe- 
nomena, this behavior can be characterized by a set of critical 
exponents of which two are independent. 

Of interest to us is the so-called percolation strength 
P{p) {p ~ Pc)^^" defined for p > pc as the probability 
that a randomly chosen site belongs to the percolating cluster, 
the average cluster size S{p) ~ \pc — p\~^^''^ , and the corre- 
lation length £^{p) ^ \pc — pl^" . We give some quantities a 
subscript "per" to avoid confusion later on, where similar, but 
different quantities appear The percolation strength is analo- 
gous to the magnetization in spin models, while the average 
cluster size corresponds to the magnetic susceptibility. 

The number density of clusters having size s, i.e., with s 
occupied sites, is assumed to be distributed according to 



ls{p) (x s ^exp(— cs), 



(1) 



where the coefficient c vanishes with an exponent 1 /a when 
the percolation threshold pc is approached from below c oc 
(Pc — pY^" ■ criticality, the cluster distribution becomes 
ls(j>c) oc s^'^ and falls off algebraically. This factor measures 
the configurational entropy of clusters, while the exponential 
is a Boltzmann factor Together with the so-called Fisher ex- 
ponent T, the exponent cr determines the critical exponents, 
such as /3per, 7pci and v through scaling relations. 

In terms of the cluster distribution 1^, the average cluster 
size is given by 



lls^^ls{p) ^ 



(2) 



where only finite clusters are included in the sum. To under- 
stand that S thus defined indeed is a measure for the average 
cluster size, note that the combination slg (p) in the denomi- 
nator is the probability that a randomly chosen site belongs to 
a cluster of size s. When summed over s, this gives the prob- 
ability that the randomly chosen site belongs to a cluster of 
arbitrary (finite) size. The ratio slg {p)/J2s (p) is therefore 
the probability that the cluster to which the randomly chosen 
site belongs is of size s. When this ratio is multiplied with the 
cluster size and summed over s, one obtains indeed a measure 
of the average cluster size. With the explicit expression for S 
one easily finds that 



an occupied site belonging to the same cluster as the origin 
does. At criticality, it has in d space dimensions the algebraic 
behavior 
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X' 



(5) 



with 77pcr a critical exponent. 

To obtain a similar expression for the percolation strength 
as given in Eq. (El) for S, one uses the identity 



p = P(p) + ^s/,(p), 



(6) 



stating that an occupied site either belongs to the percolating 
cluster or to a finite one. Hence, for p > Pc 



p{p)--J2sUp). 



and 



T-2 



(7) 



(8) 



Below the percolation threshold, P{p) = 0, so that the iden- 
tity becomes p = J2s sls{p)- 

Given the analogy with spin systems, it should be possible 
to define a partition function Z from which the percolation 
strength (magnetization) and the average cluster size (suscep- 
tibility) can be obtained by differentiation. Indeed, with the 
definition 



ln(Z) cx Is 



one finds 



dc 



(9) 



(10) 



as required by the analogy. Physically, ln(Z) denotes the total 
number density of finite clusters. 

The two critical exponents 7pcr and /3pcr can be related to 
the critical exponent a characterizing the algebraic behavior 
of ln(Z) ~ \pc — pp^" close to the percolation threshold as 
follows 



2-a 



T- 1 



2/3p 



-7pc 



(11) 



7pcr 



(3) 



For later reference, recall that the susceptibility also equals 
the integral over the correlation function G'pcr(a;) of the sys- 
tem under consideration. Because of Eq. (^ one finds for 
percolation 



^ Gper(a:) cx ^ s'^lsip), 



(4) 



where the sum J2x lattice sites. Physically, G'pcr(a;) 

is the probability of finding at a distance x from the origin 



On dimensional grounds, one expects that near the percolation 
threshold ln(Z) ^ ^^"^ ~ \p — p^'^'^ , so that one arrives at the 
hyperscaling relation 



dv = 2 



(12) 



involving the dimensionality d of the lattice. In the context of 
thermal critical phenomena, relations ( pi] ) and ( 12 ) are known 
as Rushbrooks' and Josephson's scaling law, respectively. A 
hyperscaling relation in general breaks down at the upper crit- 
ical dimension. Beyond the upper critical dimension du, the 
critical exponents are locked in their values at o? = du- For 
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percolation, as we will see below, the upper critical dimen- 
sion is = 6. 

Another scaling relation between the various critical expo- 
nents is obtained from the expressions (0) and (||) involving 
the correlation function. Using ^ as lower cutoff in the inte- 
gral, one finds 



7pcr = iy{2 ~ 7/por), 



(13) 



which in the context of thermal phase transitions is known 
Fisher's scaling law. 

The correlation length exponent v is known to have a ge- 
ometric meaning, being related to the Hausdorff, or fractal 
dimension D at criticality. The latter can be defined through 
the average squared distance between the sites of a cluster. 



(14) 



i,3 = l 



with Rs the so-called radius of gyration, and x = 
(1/s) '^he center of mass of the cluster with its sites 

at Xi, as 



Rs 



(15) 



for large enough clusters. A standard definition of the corre- 
lation length in terms of the correlation function G{x), with 
X — Ixl, reads 



.2 E.^'G(x) 



(16) 



It measures the average distance between two sites in the same 
cluster. Because of the observation (^, this definition can be 
transcribed in terms of the cluster distribution as 



Considering Eq. (Q) at criticality with a long-distance cutoff 
Xmax — L, corresponding to s^ax ~ ^ we obtain a rela- 
tion between the fractal dimension and the critical exponent 
r^pcr characterizing the algebraic behavior (||) of the correla- 
tion function. 



= i(c^ + 2-r;p„.)• 



(21) 



The de Gennes' estimate ( pO| ) is equivalent to setting r/pd to 
zero here. 

To determine the upper critical dimension it is prudent to 
consider uncorrected percolation on a Bethe lattice, which 
is exactly solvable. Because of the absence of closed paths 
on such a lattice, it mimics an ordinary lattice with high di- 
mensionality, where one expects the same critical behavior as 
in the upper critical dimension. The cluster distribution on a 
Bethe lattice is given by Eq. ^ with r = 5/2 and a = 1/2, 
corresponding to the critical exponents 

a = -1, /?per = 1, 7per =1, = 1/2, r/per = 0, 

(22) 

and the fractal dimension D ~ A. These values are consistent 
with the scaling laws in d = 6, which can thus be identified as 
the upper critical dimension. 

To summarize, percolation near criticality is specified by 
the two exponents a and t parameterizing the cluster distribu- 
tion, both of which are related to the fractal dimension D via 
Eqs. ( [l8| ) and (|l9|). Given the value of t and a, the critical ex- 
ponents can be obtained using scaling relations. Specifically, 



T - 1 

a = 2 , 13, 



riper- = 2 + d 



per 



r-2 

a 
T - 1 

da ' 



3-r 



D = 



T - 1 



(23) 



(17) 



With Eq. (p^, the relation between the correlation length and 
the fractal dimension is then easily found to be 



D 



(18) 



Combined with this relation, the hyperscaling relation ( p^ 
leads to the well-known expression for the Fisher exponent 
T in terms of the fractal dimension: 



(19) 



Using a Flory-typeD of argument known from polymer physics 
and the theory of self-avoiding walksO, de GennesB estimated 
the fractal dimension for uncorrelated site percolation to be 



DdG = (d + 2)/2. 



(20) 



In Fig. [l|, the Fisher exponent r with this estimate is compared 
with numericalid = 3,4,5) and analytic (d = 2,6) results 
taken from Ref.B. 



III. IDEAL BOSE GASES 

In this section we wish to point out a close analogy between 
percolation and Bose-Einstein condensation in an ideal Bose 
gas. The partition function Z describing this system in d space 
dimensions can be written as an integral over momentum asEl: 



HZ) 



-V 



d'^fc 

{2nY 



■In 1 



-E(k)/ksT 



(24) 



ignoring the k — Q contribution. In Eq. (|4|), V is the vol- 
ume of the system, while E{k) — k"^ /2m — pi is the single- 
particle spectrum, with m the mass of the particles and /i the 
chemical potential which accounts for a finite particle num- 
ber density. In deriving this, the time coordinate was analyt- 
ically continued to the imaginary axis, where it becomes a 
cyclic variable taking values only in the interval [0, Ti/fceT], 
with fee Boltzmann's constant. After a Taylor expansion of 
the logarithm, the partition function-takes the familiar form of 
a fugacity [~ cxp(/i//cBr)] seriestlj 



HZ) 



V 



(25) 
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distance x from the origin after w steps, while Gio(O) denotes 
the probability for the random walk to return to its starting 
point after w steps. 

As last quantity we need the correlation length ^ which is 
given by = X^k^Tj (—Anfi), i.e., ^ cx /i^^^^, as also fol- 
lows from definition ([iq). The chemical potential is negative 
in the normal state and vanishes on approaching the critical 
temperature Tc, which is determined by the condition 
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G(0)| 



(30) 
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FIG. 1: The Fisher exponent r and the exponent a parameterizing 
the cluster distribution of uncorrelated site percolation as functioHi 
of the dimensionality d of the lattice. The points represent resultsu 
obtained by lattice simulations (d — 3, 4, 5) and analytic methods 
(d = 2, 6}, while the line in the top panel is based on de Gennes' 
estimate ( pcj ) which can be given for any dimension 2 < d < 6, not 
just integer values. 



where 



IwiT) = w cxp(/^w/fcBT), 



(26) 



with T = d/2 + 1 and A = h^y2Tr/mkBT the de Broglie 
thermal wavelength. 

The correlation function G{x) written as a sum over w 
reads 



G{x) := '^Gni{x) = ^^w^^cxp ( - 



X w 



so that 



and 



J d'^xG{x)=J2' 



(27) 



(28) 



(29) 



where the left side of the last equation denotes G(0), with 
G{k) indicating the Fourier transform of G{x). Physically, 
the correlation function Gw{x) denotes the probability for a 
Brownian random walk starting at the origin to end up at a 



where n is the particle number density. This equation is valid 
for any /i < and thus for any temperature T > T^. How- 
ever, it starts to break down at the critical temperature where 
the condensate starts to form in the fc = state, which has bee 
ignored in Eq. (p4|). The critical temperature is in other words 
the lowest temperature where this equation is still valid. From 
the criticality condition it follows that upon approaching the 
critical temperature from above, the chemical potential tends 
to zero in a d-dependent way as ^(T) -(T - Tc^/^'^^'^h 
Below this temperature, in the condensed state, the chemi- 
cal potential remains zero. The d-dependence here distin- 
guishes an ideal Bose gas from a Gaussian theory, where /i 
tends to zero as ii{T) ^ \T — Td irrespective of the di- 
mensionality. The critical exponents of an ideal Bose gas, 
which can be extracted from the information given above, are 
(3 = 1/2, u = l/{d - 2), and 7 = 2/{d As noted 

a long time ago by Gunton and Buckinghamllil, these expo- 
nents are the same as for the spherical model in d dimensions. 
This model corresponds to the limit ri ^ cx3 of the 0(ri) spin 
model. 

They continued to show that an ideal Bose gas with the 
modified energy spectrum e(fc) cx is in the same uni- 
versality class as the spherical model with long-range inter- 
actions considered by JoyceU. The critical exponents for 
D < d < 2D, with d = 2D the upper critical dimension, 
are given by 



d-2D 

d-D ' 
1 



d-D' 



1 

2' 

2-D. 



D 



d-D' 



(31) 



The critical exponents of the standard ideal Bose gas are re- 
covered by setting D = 2. As will be shown shortly, the 
energy spectrum index D equals the fractal dimension. 

The partition function of the ideal Bose gas with a modified 
energy spectrum can again be represented in the form (|2^), 
with T = d/D + 1 and fi (xjT - Ty^", where a = l/vD 
as for percolation [see Eqs. (jlq) and (tt9h]. The thermal wave- 
length A for the generalized Bose gas is defined in a way that 
the spectrum can be written as 



e/ksT = coX^'k 
with the constant cd conveniently chosen such that 



d/D 
=15 



1 



2 T{d/D) 



(47r)d/2 D r{d/2) 



(32) 



(33) 
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FIG. 2: A worldline loop involving three particles projected onto the 
finite imaginary time interval < r < h/k^T (left panel). The 
loop, which has winding number u; = 3, can be equivalently repre- 
sented as the worldijnes of three cyclically permuted particles (right 
panel). (After RefB.) 



where T{x) is the Euler gamma function. The correlation 
function G{x) at x is again given by Eq. ( |2^ while also 
the integral of this function over space yields the same result 
as for D ~ 2 given in Eq. (p9[). From the criticality condition 
( |3C| ) it follows that the correlation length scales with the chem- 
ical potential as ^ cx /i^^^^, where — as already indicated by 
our choice of notation — D denotes the fractal dimension as 
calculated from its definition (|l5]). The critical exponents (|I]) 
now follow easily. _ 

As detailed in Ref.lHl, the fugacity series (^Sj) represents a 
sum over loops with the distribution (Eq), or, equivalently. 



(34) 



with £ = wa the "length" measured in units of the charac- 
teristic "length" scale a = h0, and 9 = —^k-oT/h the string 
"tension". The first factor in ( |34| ) measures the configurational 
entropy of loops, while the second is a Boltzmann factor. Each 
loop is characterized by the winding number w, telling how 
often it wraps around the imaginary time axis. Physically, a 
loop represents the worldlines in imaginary time of w particles 
grouped together in a single ring known from Feynman's the- 
ory of the A-transition in ''Heo. Particles in the same ring are 
cyclically permuted after an imaginary time h/k^T (see Fig. 
g for an illustration involving three particles). With increas- 
ing imaginary time, the particles move along closed strings 
towards the initial position of the particle in front of them (see 
Fig. U). When only one particle is contained in a ring, the 
particle returns to its own initial position after an imaginary 
time h/k-QT. On approaching the critical temperature from 
above, the chemical potential and thus the string tension be- 
come smaller, and loops with larger winding numbers start to 
appear. At Tc, the string tension vanishes and long loops con- 
taining arbitrarily many cyclically permuted particles appear 
in the system, signalling the onset of Bose-Einstein condensa- 
tion. Above the critical temperature, long loops are exponen- 
tially suppressed. Because the loops are worldlines embedded 




FIG. 3: Worldline loops of 20 particles in two space dimensions. 
Some of the worldlines are grouped together in single rings. The 
beats in the top panel represent particles moving in imaginary time 
along the strings towards their own initial position, in case of a single 
particle in a ring, or towards that of the neighboring particle in front 
of them. After an imaginary time h/ksT, these positions are reached 
and the particles cyclically permuted. In the bottom panel the beads 
are omitted, leaving behind a loop gas of worldlines. (After Ref.u) 



in spacetime and parameterized by Euclidean time, £ and a 
have the dimension of time, not of length — whence the quota- 
tion marks used below Eq. (^4|). 

From this and also from comparing the partition function 
( |2^ ) of an ideal Bose gas with the one for percolation, Eq. (^), 
we observe a close connection between the two phenomena. 
The rings of cyclically permuted particles in a Bose gas cor- 
respond to the clusters of percolation theory. The number of 
Bose particles contained in a ring corresponds to the number 
of sites contained in a cluster. In the same way that clusters 
grow when the percolation threshold is approached from be- 
low, rings with larger winding numbers appear in a Bose gas 
on approaching the critical point from above. In both cases, 
the behavior at the critical point is algebraic, rather than expo- 
nential. The percolating cluster atp>pc corresponds to the 
Bose-Einstein condensate at T < Tc. [The reason for ignor- 
ing the k = contribution in Eq. ( p4| ) was to bring out this 
analogy as clearly as possible.] But there are also differences. 

Although the exponents ( plf ) satisfy the scaling relations, 
(3, 7, and 77 are not related to the exponents cr and r parame- 
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terizing the loop distribution ( [26[ ) in the same way as in per- 
colation theory [see Eq. (^3|)]. Using these equations with the 
ideal Bose gas values for r and a, we find instead 



Ppcr — 1; Tpc 



d-2D 
' d-D ' 



Tjp,, = 2 + d-2D, (35) 



where we have given the exponents also the subscript "per" to 
indicate that we used the percolation formulas for them. 

Despite the differences in form, these exponents too satisfy 
the scaling relations 2/3por + 7Dcr = 2 — a and 7por — v{2 — 
Vpcr), with a and v given in (|T]). Apparently, there are two 
different sets of critical exponents, both satisfying the scaling 
laws. 

To understand this, let us compare the sum rules and 
( |2^ ) satisfied by the correlation function of percolation and an 
ideal Bose gas, respectively. There is no obvious connection 
between the two, implying that the "susceptibility" and also 
the "magnetization" are defined differently in the two systems. 
In percolation theory, they were defined in terms of the parti- 
tion function in Eq. (p^. With c replaced by — /i/fceT to map 
the cluster distribution of percolation theory onto the loop dis- 
tribution of an ideal Bose gas, it follows that P and S corre- 
spond to the condensate density no and the compressibility k 
of an ideal Bose gas, respectively, with 



no 



(Tc - T) 



/3pc 



\Tr. - T\ 



(36) 



The exponents f3 and 7 featuring in Eq. ( |3l| ) are, on the 
other hand, connected to the order parameter of an ideal Bose 
gas, which is the field (po describing the condensate, with 
I '/'op = '^0- This last equation relates the two /9-exponents, 
giving /3por = 2(3. Since /9per = 1, it also follows that 
7pcr — ^Oi, as one would expect naively. It should be noted 
that the exponents v and a are not sensitive to the nature of 
the order parameter describing the transition. Given the pres- 
ence of a single diverging length scale, any sensible definition 
of the correlation length leads to the same exponent i'. And 
a merely characterizes the algebraic behavior of the logarithm 
of the partition function, which is essentially the pressure. De- 
spite the difference in the correlation function of percolation 
and an ideal Bose gas, the relation ( p^ ) between the correla- 
tion length exponent i' and the fractal dimension is the same 
for the two, as can be easily checked repeating the derivation 
of this relation with the correlation function of an ideal Bose 
gas. 

A more fundamental difference between percolation theory 
and the string picture of an ideal Bose gas is that in the latter, 
the number of particles is fixed, while in the former there is 
no condition on the number of occupied sites. The fixed num- 
ber of particles in an ideal Bose led to the criticality condition 
( |3C| ) with G(0) given in Eq. ( |2^ which in turn determined the 
critical exponent ly. In percolation theory, where there is no 
such a condition, ly is not determined that way. As a result, 
although the two phenomena are described by essentially the 
same theory specified by a loop or cluster distribution, they in 
general are not in the same universality class. Only the spe- 
cial case of percolation in the upper critical dimension d^ — 6, 
where Z) = 4, or on a Bethe lattice is in the same universality 



class as an ideal 6-dimensional Bose gas with fractal dimen- 
sion D = A. 

However, in the absence of a criticality condition, we ex- 
pect uncorrelated percolation and strings of the type appearing 
in the description of an ideal Bose gas to be described by ex- 
actly the same theory. This is indeed what has beenjObserved 
in numerical simulations on random string networkstJ. 



IV. RANDOM STRING NETWORKS 

These simulations describe theories with a global U(l) 
gauge symmetry, which undergo a phase transition from a 
symmetric state to one where the U(l) symmetry is sponta- 
neously broken. The ordered state, which is characterized by a 
complex order parameter, can have topological line defects — 
known as cosmic strings in the context of cosmology, and as 
vortices in the context of condensed matter. The line defects 
are either closed or terminate at the boundary, they cannot 
terminate inside the system. Upon circling a vortex of unit 
strength, the phase d of the order parameter changes continu- 
ously by 2tt. 

To npperically study line defects, Vachaspati and 
Vilenkiro considered a cubic lattice with periodic boundary 
conditions. They discretized the vacuum manifold of the or- 
dered state by allowing -d to take on only certain discrete val- 
ues in the interval < d < 2tt. The restriction to this interval 
assures that only vortices of unit strength are generated. Each 
lattice site is assigned one of the discrete values at random. 
Phases at different lattice sites are therefore uncorrelated. Vor- 
tices are traced by going around each plaquette of the lattice 
once. If a phase difference of 2tt is found upon returning to 
the site one started at, it is concluded that a vortex penetrates 
the plaquette. In going from one lattice siteJo a neighboring 
one, the geodesic rule introduced by KibbleO is implemented, 
assuming that the (discrete) changes in d are minimal. When 
two instead of just one vortex is found to penetrate a unit cell, 
the two incoming and two outcoming vortex segments are ran- 
domly connected. There are different ways to do this, the sim- 
plesj-one being to choose equal probabilities as was done in 
Ref.E3. But there are other choices, such as favoring those 
possibilities that generate longer loopsEJ. In this way a whole 
random network of vortices (of unit strength) can be traced 
out. 

Such a random string network is of relevance to the forma- 
tion of topological defects after a rapid phase transitipn in- 
volving the spontaneous breakdown of some symmetryllj. As 
was first proposed by KibbleO in the context of cosmology, 
causally disconnected pieces of the early universe may after a 
rapid temperature quench end up in different, but degenerate 
states of the broken-symmetry vacuum, characterized by dif- 
ferent values for the phase of the order parameter. When 
these regions come together, the resulting frustration at the 
boun 
Zurel 



ries may lead to the formation of topological defects. 



noted that a similar mechanism might operate in the 
context of condensed-matter physics and developed a qualita- 
tive theory, based on scaling arguments, of the formation of 
topological defects in time-dependent phase transitions. To 
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test this scenario, superfluid ^He-B in rotation was radiated 
with neutronsEl. After a neutron-absorption event, in which a 
small region is briefly heated up to the normal phase, a vortex 
network was produced. Under the action of the Magnus force, 
some of the vortex loops in the network expand and connect 
to the container wall. They are then, still under flie action of 
the Magnus force, pulled straight and, finally, end up in the 
center of the cylindrical container aligned parallel to the rota- 
tion axis, where they can be detected by means of NMR. The 
number of vortices produced in the neutron irradiation exper- 
iments was in agreement with Zurek'S|-flredictions. The vor- 
tex network was numerically simulatedcl with the Vachaspati- 
Vilenkin algorithm discussed above. 

Besides finite vortex loops, also infinite vortices appear in a 
random network with high enough vortex densities. The divi- 
sion between finite vortex loops and infinite ones is of course 
somewhat arbitrary on a finite lattice. A possible choice is to 
classify vortices with a length larger than 2L^, where L is the 
lattice size, as infinite. Universal quantities such as the critical 
exponents should be independent of the precise jdnoice. The 
fractal dimension found in the initial simulationslij were con- 
sistent with that of a Brownian random walk, D — 2. This 
result can be understood^ by realizing that a dense random 
network was simulated, where a single vortex experiences an 
effective repulsion from the neighboring vortex segments as 
the volume occupied by these segments is not available any- 
more. The situation is similar to a polymer in a dense solution, 
which also displays the structure of a Brownian random walk. 

Without any external field, all phase values of the order pa- 
rameter have equal probability. VachaspatiH lifted this degen- 
eracy of the vacuum manifold by giving some values a bias, 
thus obstructing the formation of vortices and decreasing the 
vortex density. When studying the statistics of the random 
network as function of the bias, a threshold was found, above 
which infinite vortices are absent and finite loops exponen- 
tially suppressed. The sudden appearance of infinite vortices 
as the threshold is approached from the low-density side, was 
interpreted as a percolatinn-process. 

Strobl and HindmarshEj carried out more extensive nu- 
merical simulations on the transition of a random string net- 
work triggered by changing the bias. Using the percolation 
set of critical exponents, they measured — to within statisti- 
cal errors — values similar to those for 3-dimensional uncor- 
rected site percolation, which has cr = 0.45, t = 2.18 and 
D = 2.53Q. For the vortex loops, the value D > 2 means that 
they are self-seeking. The conclusion that the phase transition 
of a random string network is identical to the uncorrected per- 
colation transition is precisely what one would expect given 
our observations of the last section. 



V. CORRELATED PERCOLATION 

Up to this point we considered only uncorrelated site per- 
colation, where lattice sites were occupied at random, and 
nearest-neighboring occupied sites grouped together in clusp 
ters. However, as is known from work on the Ising modelQ, 
this is not the whole story. 



In the context of the Ising model, it is more appropriate to 
refer to occupied sites as sites with spin up, say. In general, 
the percolation threshold temperature is below the critical one. 
That is to say, already before the thermal phase transition is 
reached, clusters of spin-up sites percolate the lattice. Only 
on a 2-dimensional square lattice, the site percolation thresh- 
old temperature coincides with the critical temperature. De- 
spite this, the Ising-model and percolation critical exponents 
are different-as in higher dimensions. The situation is similar 
to the one discussed at the end of Sec. |l|, with the percolation 
transition being in a different universality class as the thermal 
one. 

It is nevertheless possible to modify the definition of clus- 
ters in such a way that the critical behavior of the modi- 
fied clusters becomes identical to that of the Ising modelEj. 
This apprpach goes back to seminal work by Fortuin and 
KasteleyncJ, who mapped the partition function of the Ising 
model onto a correlated percolation problem. To obtain a 
modified cluster, take two nearest-neighboring sites of 
a standard cluster with all the spins up, and add the pair 
to the new cluster with the bond probability pij = 1 — 
exp(— 2 J/fceT), where J is the spin-spin coupling of the 
Ising model. The factor 2 J appearing here corresponds to 
the increase in energy when one of the two spins involved in 
the bond is flipped. The upshot is that the bond clusters are in 
general smaller than the standard ones and also more loosely 
connected. The bond clusters percolate right at the critical 
temperature and display critical behavior identical to that of 
the Ising model. 

The method has been turned intcua powerful Moate Carlo 
algorithm by Swendson and WangE3, and by WolflEj, where 
not individual spins are updated, but entire clusters. The main 
advantage of the nonlocal cluster algorithm in comparison to 
a local update algorithm, like Metropolis or heat bath, is that it 
drastically reduces the critical slowing down near the critical 
point. 

Very recently, also the critical behavior of the 3- 
dimensional 0(n) spin models, with n =,-2,3,4, was de- 
scribed as a correlated percolation procesfJ. The relevant 
clusters percolate at the critical temperature, and were shown 
to have the same critical exponents as the corresponding mod- 
els. 



VI. WORLDLINE LOOP GASES 

In the same way that correlated percolation can describe 
thermal phase transitions of interacting systems also corre- 
lated worldline loop gases can. Recall that the loop gas de- 
scription of an ideal Bose gas hinged on the representation 
( 25 ) of the partition function in terms of the loop distribution 
(34), which derived from the representation ( |24| ) as an integral 
over momentum. If we assume that near the critical point, the 
partition function of an interacting system can be put in a sim- 
ilar form, but with E{k) replaced by the single-particle spec- 
trum of the excitation that becomes critical at T^, the integral 
can in principle be done, yielding some loop distribution 
with w physically denoting the number of times a worldline 
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loop wraps around the imaginary time axis. As for the ideal 
Bose gas, this distribution is parametrized by two exponents 
a and t, and can be put in the general form (^4[), with values 
characteristic for the universality class under consideration. 
The string tension of the worldline loops vanishes at Tc, and 
loops with arbitrarily large winding number w appear in the 
system, signalling the onset of Bose-Einstein condensation in 
the low-temperature phase. Note that the representation ( p4| ) 
of the partition function as an integral over momentum does 
not require a nonrelativistic theory. Also relativistic theories 
lead to a similar fortH, for the (equilibrium) partition function 
at finite temperatureEa, as it is essentially fixed by Bose statis- 
tics. 

The generalization of the sum rule to correlated world- 
line loop gases reads 



-0l/ksT 



(37) 



where we recall that the string tension behaves near the critical 
point as 61 cx (T - T^Yl". The sum rule leads directly to the 
relation 



7 = l/o-, 



(38) 



while at criticality, with the long distance cutoff Xmax — L, 
corresponding to tUmax — as before, it gives the relation 



2-D. 



(39) 



These two relations generalize the corresponding results in 
Eq. ( pH ) for an ideal Bose gas with modified energy spectrum. 
The value 77 = separates self-avoiding worldlines {D < 2) 
from self-seeking ones. Since the fractal dimension has to be 
positiiie and smaller than the dimension d of the embedding 
spaceEj, 2 ~ d < rj < 2. Translated in terms of limits on 
T, it follows that 2 < r < d/2 + \ for self-seeking world- 
lines, while T > d/2 + 1 for self-avoiding ones. Note that in 
two space dimensions and below, no self-seeking strings are 
possible. 

As in Sec. the exponents here are those for describing 
the thermal phase transition of a loop gas of strings (world- 
lines). They are to be distinguished from the ones defined in 
the context of percolation theory given in (^3]). The two sets 
have, however, the exponents a and v in common, implying 
that the present set can also be expressed entirely in terms of 
the exponents a and r parameterizing the loop distribution 
0. Indeed, we find 

„ T- 1 ^ T-2 1 

« = 2 , /3=^^, 7=-, 

r - 1 d 
v=— — , D= -. (40) 



77-2- 



T-V 



da 



T - 1 



Since only general scaling laws are used in deriving these re- 
sults, we believe they hold for any critical theory specified by 
a loop distribution ( |34| ) with exponents a and r. It can be eas- 
ily verified that for an ideal Bose gas with a = d/ D — 1 and 
T = d/ D + \, Eq. (ES) reduces to the correct values (pT|). 



Random walk 


(7 r 


a j3 "f f) V D 


Brownian 
Smooth 


1 l + d/2 
1 l + d 


2 - d/2 d/4 - 1/2 1 1/2 2 
2-d d/2 - 1/2 1 1 1 1 



TABLE I: Critical exponents for Brownian and smooth random 
walks, together with the exponents a and r parameterizing the loop 
distribution (|34|). 



As an illustration, let us discuss the worldline loop distri- 
bution of some well-known statistical models. We start with 
random walks. In Table |, the critical exponents of two dif- 
ferent types of random walks are given, together with the cor- 
responding values for a and t parameterizing the loop distri- 
bution (|34|). Both have the value a in common. The 
first one is the Brownian random walkCl, which has a frac- 
tal dimension D = 2_isihile the second one is the so-called 
smooth random walkoO, which has D — 1 and is therefore 
self-avoiding. Whereas a Brownian random walk can be un- 
derstood as a particle hopping through a lattice with its next 
(nearest-neighboring) site chosen randomly at each step, a 
smooth random walk can be understood as a particle hopping 
through the lattice with its velocity rather than its position be- 
ing changed randomly at each step. As a result, whereas a 
typical path of a Brownian random walk is continuous but not 
differentiable, a typical path of a smooth random walk is dif- 
ferentiable while only the first derivative is not. The reluc- 
tance to change directions is typical for fermionic particles. 
To understand the value rj — 1 for a smooth random walk, 
note that the correlation function of a Brownian random walk 
is given in the continuum by. 



Gix) 



{2nY 



with 



GB(fc)= 



(41) 



(42) 



leading to the algebraic behavior (||) at criticality with rj = 0. 
A smooth random walk, on the other hand, haSpa-jGorrelation 
function involving one factor of momentum lessEaLlI, 



Gs(fc) = 



(43) 



where e is some unit vector At criticality, this again gives the 
algebraic behavior (jsj), but now with rj — 1. 

In two space dimensions, the critical exponents of various 
models are known exactly. For example, for 0(n) spin mpidels 
with n = —2 cos(27r/t) and I < t < 2 they are given byH, 



2t-3 
1~2' 



2(i-2)' 



(^-l)(^-3) 
%t{t - 2) ' 

2t 
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D = 



e + 3 

'4t(i- 2)' 



2t 



(44) 
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Model 


n 


a T 


a /3 7 T) V D 


Gaussian 


-2 


1 2 


10 1 1/2 2 


SAW 





32/43 91/43 


1/2 5/64 43/32 5/24 3/4 43/24 


Ising 


1 


4/7 15/7 


1/8 7/4 1/4 1 7/4 


XY 


2 


15/7 


1/4 7/4 


Spherical 


oo 


2 


1/2 2 



TABLE II: Critical exponents for the 2-dimensional 0{n) spin 
model, with n = —2, 0, 1, 2, oo, respectively, together with the ex- 
ponents a and r parameterizing the loop distribution (M). 



leading to D' = 7/4 indeed fulfill the relations (|6|). 

Recent high-precision Monte Carlo, studies of the 3- 
dimensional XY-universality class gaveEj 

1/ = 0.6723(3), 77 = 0.0381(2). (48) 

The resulting exponents parameterizing the worldline loop 
distribution are 

cr = 0.7582, r = 2.5291, (49) 

corresponding to a fractal dimension D = 1.9619 very close 
to two, in accordance with the smallness of -q. Since D < 2, 
the worldline loops of the 3-dimensional XY-model are self- 
avoiding. 



corresponding to 



VII. VORTEX LOOP GASES 



a = -'-^, r^(l±lM±^. (45) 



t2 + 3 



t2 + 3 



For special values of the parameter t such that n is an integer, 
one obtains the critical exponents for the Gaussian model, or 
the Brownian random walk (n = —2), the self-avoiding walk 
(SAW) (n 0), the Ising model (n = 1), and the XY-model 
(n = 2). The critical exponents for these models together with 
those of the spherical model are collected in Table O. 

A few points may be worth noting, (i) Both the Ising 
(t = 3/2) and the XY-model {t = 2) have the same fractal 
dimension. This is made possible by a minimum in D{t) at 
t = \/3 located between the two models, where D = 

(ii) In the literature, the fractal dimension is often equated 
to 1/i^ rather than to 1/aiy as in Eq. (jl8|). We believe, the 
result (|l8|) is more general. 

(iii) The special character of the 2-dimensional XY and 
also of the spherical model is seen in the loop gas picture by 
the vanishing of the exponent a, implying the absence of the 
Boltzmann factor in the loop distribution function (Q). 

(iv) Both the 2-dimensional Ising and XY-model share 
the same r-value and their worldlines have, consequently, 
the same fractal dimension. There is at least one other 1- 
dimensional object embedded in two space dimensions with 
the same fractal dimension. Surprisingly, it is the perimeter 
of a percolation cluster of standard, uncoiTelated percolation 
in two dimensionsLJ, which can also be mapped onto various 
types of random walks. Close to critic ality, the perimeter of 
the percolating cluster is again distributed according to Eq. 
(|l|), but with different values for the exponents. To avoid con- 
fusing with the 2-dimensional cluster exponents a = 36/91 
and T = 187/91 corresponding to D ^ 91/48, we call the 
perimeter exponents a' and r'. Since there is only one correla- 
tion leugth in the system, the cluster and perimeter exponents 
satisfya 

(r- 1)/ct = (r' - l)/cr', a/a'^D'/D. (46) 
The perimeter valuesS 

a' = 3/7, t' = 15/7, (47) 



Let us next turn to interacting systems with a global U(l) 
gauge symmetry, which have vortices as topological defects 
in the ordered state. Thermal phase transitions in systems of 
this type have, besides the conventional Landau description, 
an alternative description in terms of proliferating^yorticesLJ. 
Such a scenmo was first put forward by Onsagenil and later 
by Feynmanta to describe the A-transition in ^He. In the or- 
dered, superfluid state, where the U(l) symmetry is sponta- 
neously broken, only finite vortex loops are present. Because 
of the finite string tension, the loops are exponentially sup- 
pressed. Upon approaching the critical temperature from be- 
low, the string tension vanishes, and vortices start to prolifer- 
ate. For vanishing string tension, vortices can gain configu- 
rational entropy by growing without energy cost. The result- 
ing infinite vortices disorder the system and thereby restore 
the spontaneously broken U( 1 ) symmetry. The importance of 
vortex loop excitations in triggering the phase transition has 
since then been emphasized by various authorsE3. Early nu- 
merical evidence for this picture based on Monte Carlo simu- 
lations of thp-jJ-jdimensional XY-model was given by Janke 
and KleinertEJo. Analytic methods to describe the transi- 
tion using vpitex loops in a qualitative way were developed 
by Williamstll, and by Shenoy and collaborators^. 

Like worldline loops, vortex loops are specified by a loop 
distributioH, as in Eq. (|^ parameterized by two exponents 
cr' and t'cJ. These determine the critical exponents of the 
system again through Eq. (|40|), with a and r replaced by a' 
and t'. The result that the exponent l/cr' parameterizing the 
vanishing of the vortex string tension directly gives the critical 
exponent 7 — and not 2v ^^if often claimed in the literature. 



also by thi 
and Sudbc 



□resent author 



-was first observed by Nguyen 



while the connection rj = 2 — D foiLvortex loops 
was very recently derived independently in Ref.E3. 

In the past few years various groups have continued 
to numerically investigated the 3-dimensional XY-model 
from the nejjs«ective of vortices, in particular their loop 
distributionll3c3c3. Let us discuss the numerical values for the 
exponents cr' and r' obtained in these studies from the present 
perspective, concentrating on insights not available before. 



10 



Irrespective of the description chosen, either with the help 
of worldline loops or vortex loops, the expression for the cor- 
relation length exponent is formally the same. Since also the 
numerical value of this critical exponent is independent of the 
description, the exponents a and t parameterizing the world- 
line loop distribution and a' and r' parameterizing the vortex 
l oop distribution satisfy (r — l)/^ ~ (r' — 1)/(t' as in Eq. 
(pq). They should both lead to the same numerical value fot-j/ 
given in Eq. (Bsh. Using for a' and r' the numerical valuesLJ 



1/cr' = 1.45(5), t' = 2.4(1) 

300, 



(50) 



reported by Nguyen and Sudb0EI, we find from Eg. ( [40[ ) a 
result V = 0.68 very close to the expected one (Q). Since 
2 < t' < d/2 + 1, vortex loops in the XY-model are (slightly) 
self-seeking. The numerical values for a' jatid r' reported |b* 
the other groups give the results v ~ 0.62CJ and v = O.60IIJ, 
which are within about 10/Lpf the expected result. 

It was noted by WilliamsEl that the fractal dimensioruD ex- 
tracted from the value r' — 2.23(4) reported in Ref.a was 
close to the estiniate D = {d + 2)/2 obtained from a Flory- 
type of argumenlH. This estimate is identical to de Gennes' 
one (20) for uncorrected percolation which was also based 
on a Flory-type argument. It should be noted, however, that 
de Gennes' modification of the original Flory argument was 
motivated in part to arrive at the correct upper critical dimen- 
sion ~ 6 for percolation, rather than = 4 for random 
walks and the XY-model. |_. 

Kajantie and coUaboratorsllj investigated different algo- 
rithms to construct the string network. Depending on the al- 
gorithm chosen, they obtained (slightly) differing quantitative 
results. None of the algorithms used yielded a percolation 
threshold precisely at the critical temperature. From this ob- 
servation they concluded that "geometrically defined percola- 
tion observables need not display universal properties related 
to the critical behavior of the system". It has been suggestedH 
that this discrepancy might instead hint at the existence of 
a unique, yet to this date unknown, algorithm to construct 
the string network that would give precisely the expected XY 
critical behavior The situation would then be similar to that 
found for correlated percolation and thermal phase transitions 
in 0(n) spin models discussed in Sec. 0. 



VIII. CONCLUSIONS 

In this paper we discussed the close analogy between per- 
colation of (correlated) clusters and proliferation of strings at 
thermal phase transitions of interacting systems. Two types 
of strings were considered: worldlines and topological line 



defects. Both are described in a way similar to clusters in per- 
colation theory, being specified close to the critical point by a 
distribution parametrized by only two exponents. At the crit- 
ical temperature, the string tension vanishes and loops prolif- 
erate in the same way as clusters percolate. If the loops are 
finite-temperature worldlines, their proliferation signals the 
formation of a Bose-Einstein condensate which spontaneously 
breaks the global symmetry characteristic of the phase transi- 
tion under consideration. If, on the other hand, the loops are 
vortices, their proliferation signals the disordering of the or- 
dered state in which these vortices are topologically stable. 

In deriving our results we used the ideal Bose gas with mod- 
ified energy spectrum as a stepping stone. The model, despite 
being noninteracting, has nontrivial critical exponents which 
are known exactly. It also has the virtue that it can be mapped 
onto a loop gas in an exact way. The resulting relations be- 
tween the exponents parameterizing the loop distribution and 
the critical exponents describing the thermal phase transition 
of the system were shown to be easily generalized to interact- 
ing loop gases using general scaling laws. As for percolation, 
the two exponents parameterizing the distribution were shown 
to encode the entire set of critical exponents. 

Various statistical models with exactly known critical expo- 
nents were discussed from the perspective of a worldline loop 
gas. Also recent numerical studies on the statistical properties 
of the vortex loop gas in the 3-dimensional XY-model were 
discussed. 

It was shown that a random string network is to thermal 
phase transitions involving correlated strings as what uncor- 
rected site percolation is to thermal phase transitions in 0(n) 
spin models. 

A possible connection of the approach discussed in this 
paper wilh the Monte Carlo loop algorithm introduced by 
Everetza is presently under investigation. 
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